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In order to analyse the role of the quasiparticle-phonon interaction in the origin 
of nuclear gap, we applied an approach which is similar to the Eliashberg theory 
for usual superconductors. We obtained that the averaged contribution of the 
quasiparticle-phonon mechanism to the observed value of the pairing gap for 120 Sn 
is 26% and the BCS-type mechanism gives 74% . Thus, pairing is of a mixed 
nature at least in semi-magic nuclei — it is due to the quasiparticle-phonon and 
BCS mechanisms, the first one being mainly a surface mechanism and the second 
one mainly a volume mechanism. The calculations of the strength distribution for 
the odd-mass nuclei 119 Sn and 121 Sn have shown that the quasiparticle-phonon 
mechanism mainly improves the description of the observed spectroscopic factors in 
these nuclei. For the case of nuclei with pairing in both proton and neutron systems 
it is necessary to go beyond the Eliashberg-Migdal approximations and include the 
vertex correction graphs in addition to the rainbow ones. The estimations for 
spectroscopic factors performed within a three-level model have shown that the 
contribution of the vertex correction graphs was rather noticeable. 



1 Introduction 

As it is known in the microscopic theory of ordinary superconductors, the BCS 
model of superconductivity is the limit case g 2 < < 1 of the Eliashberg theory in 
which the electron interaction is determined only by the quasiparticle-phonon 
coupling gEl. In the microscopic theory of finite nuclei, the BCS equation 
for the gap is rU§ed , as a rule, which contains a phenomenologically chosen 
pp interaction an. It would be justified if g 2 << 1 for real nuclei where g is 
the correspondingly dimensionless phonon creation amplitude. However, we 
have only g 2 < 1 in magic □ and in semi-magic □ nuclei and the case of strong 
coupling g 2 > 1 for nuclei with pairing in both nucleon systems. Therefore , 
one can expect a contribution of the quasiparticle-phonon interaction (QPI) 
to the observed nuclear gap. 

This question has another important aspect. As it is known, the most 
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collective low-lying phonons, which make the largest contribution to the QPI 
effects in nuclei, arc mainly surface excitations. Thus the explicit singling out 
of the quasiparticle-phonon mechanism of nuclear pairing will make it possible 
to answer the old question whether nuclear pairing is a volume or surface effect. 

And last but not the least , a consistent study of the QPI contribution to 
pairing should improve the description of nuclear excitations. At present, this 
is especially interesting in connection with the quick development and using the 
qualitatively new gamma-ray arrays like EUROBALL and GAMMASPHERE. 

The problem of the origin of nuclear pairing has been long discussed on 
a phenomenological level within the theory of finite Fermi systems using its 
idea of the internal (volume) and external (surface) pp interaction (see Q) . 
Microscopically , the question about a^contribution of the QPI to the gap 
has been discussed in 13 and recently inSoE]. InEl the authors solved a BCS 
equation with the pp interactionxaused by an exchange of a collective phonon 
between two quasiparticles. In Em the authors also considered the "insertion" 
graphs which correct single-particle moving. However, for the nuclear case of 
g 2 > 1 it is necessary to go beyond the Eliashberg-Migdal approximation which 
complicates the problem strongly. 

In this article we calculate the contribution of the QPI to the observed gap 
within the approximation of g 2 < 1 for the nucleus 120 Sn taking into account 
both types of graphs and also consider the case of strong coupling g 2 > 1 using 
a simple 3-level model. 

2 General relations 

The initial equations of our approach are very general equations for the 
one-particle Green functions in a Fermi system with pairing G, G h , F^ and 
F^ 2 ) which contain general mass operators S, Y, h , E^ and E^ 2 ). In order to 
single out the well-known components, i.e. the mean field and pairing, we 
represent each of the mass operators as a sum of two terms the first one being 
energy-independent and the second one energy-dependent: 

E(e)=t + M(e), E«(e) = P+M«( £ ), (1) 
Z^ 2 \e) = t^+M^ie) = A^+M^ 2 \e), 

where the quantities E, correspond to a mean field and E^, E< 2 > describe 
a pairing of the BCS type. The quantities M, M {2 \ contain the 

QPI and do not fix so far. 

With taking into account Eqs. (0) the general system for the Green func- 

2 



tions can be transformed into the following equations (see IS'tEl for derivation): 
G = G + GMG - F w M (h) F {2) - GM {1) F {2) - Af (2) G (2) 

p (2) = p(2) + p{2) M Q + Q(h) M (h) F (2) _ F (2) M {X) F (2) + g(/») M (2) G] 

(and the same for G w and F^). The bare Green functions G, G (h) and 
are the well-known Green functions of Gorkov which contain single- 
particle energies s\ and gap A A . 

In what follows the input parameters of our problem are phenomenological 
(observed) single-particle energies e\ and gap values A\ The latter are usually 
taken from the experiment or obtained from solving the BSC equation with a 
phenomenologically chosen pp interactionB'u. Because in the phenomenological 
e\, A\ there is a contribution of the corresponding MW terms, the quantities 
ex, A\ should be "refined" from these contributions to avoid a double counting 
of the MS l > terms or ,in our case, of the QPI. In other words , a refined basis 
i\,A\ should correspond to the mass operators and should be used in 
the calculations. The procedure. to obtain a connection between S\,A\ and 
e>, A\ has been described inBlMJ. The final formulae are the following 



i\ + Me 



l + qx(E x ) 



(3) 



= a1 1i2) = 



A A + M^' 2) (^ A ) 



where E\ = e\ + A\, q\ = -M od d\{E\) / E\ and M eveni M od d are the even 
and odd terms of the mass operator M. In equations (0) the energies i\ and 
S\ are reckoned from the corresponding chemical potential fi and fi. These 
results were obtained by using the approximation which is diagonal in the 
single-particle index A. 

Thus, in order to calculate physical characteristics it is necessary to solve 
Eqs. @ to obtain the new basis i\, A\. However , from the point of view of 
the nature of nuclear pairing the question arises about a possible contribution 
of the QPI to the quantity A\ because the difference between A\ and A\ is 
only due to the explicit singling out of the QPI contained in the terms. 
Because A\ is energy-independent one can write a general equation for it (see 
alsoffiS) 

= E W xxx^Sh ( 4 ) 

A' 
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where W is a new energy-independent pp interaction and F^ 2 ) satisfies Eq.(||) 
in which the quantities MS*' may be taken in any approximation we need. 
Further we will consider our g 2 approximation (g 2 < 1) and a model case of 
the strong coupling g 2 > f . 

3 g 2 approximation (g 2 < f) 

In this approximation we should take 



M = -Q ^_ (5) 

where the circle is the phonon creation amplitude g and the single line means 
the Green function in its matrix form which includes G , G h , p( 2 \ Here, 
the difference between A\ and A\ is due to the exchange of one phonon be- 
tween two quasiparticles (see Eq.(||)) but, in accordance with the Eliashberg 
theoryEl, this difference is corrected by the (f -I- qx)" 1 factor. 

The system @ has been solved for the semi-magic nucleus 120 Sn. At first 
phcnomcnological E\ and were obtained from the existing experimental 
data for the neighbouring 119 Sn and 121 Sn nuclei. The BCS equation for Aa 
was solved with the phenomenological volume pp interaction taken fromlia. 
The phonons have been calculated within the theory of finite Fermi systems 
H.We used 21 of the most collective 2 + , 3~, 4 + , 5~ ancL6 + phonons with the 
energy not exceeding the neutron binding energy. SeeH, where values of Ax 
and A\ are given, for detailed discussions. Here, for simplicity we use a simple 
averaging procedure, for example for A^: 

For the quantity (A av — A av )/A av , which gives the QPI contribution caused 
only by the phonon exchange between two quasiparticles (or by the retarded 
pp interaction), we obtained the value of 31%. 

Further, we have performed the calculations in the g 2 approximation for 
Ax according to Eq.(|j), i.e. we have solved this equatioii_with the Green 
function F^ 2 ) obtained from Eq.(||) in our g 2 approximation E3. The very first 
term of Eq.(|) is so that the equation A(BCS) = WF® is a BCS-like 
equation which determines a "pure" BCS part of the phenomenological gap 
value if the new pp interaction W is known. Ear simplicity the interaction 
W was taken in the same functional form as inli-J, but the parameter c p was 
determined from the condition that the average value A a „ found by the solving 
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Table 1: Spectroscopic factors in 119 Sn (first line) and 121 Sn (second line) calculated 
with (third column) and without (fourth column) taking into account the effect of the new 
(refined) basis. 



A 


n S x 


Exp.B 


R+ 


R- 


lgl/2 


0.75;0.6 


0.66(0.0-3.0) 


0.54 








2d3/2 


0.4;0.45 


0.4 


0.28 


0.44;0.65 


0.35 


0.42 


3sl/2 


0.26;0.32 


0.36 


0.23 


0.3 


0.43 


0.49 


1M1/2 


0.29 


0.26 


0.19 


0.49 


0.43 


0.45 



of the equation obtained and of the system (|3|) are the same. The results 
of these calculations are as follows. We obtained that the contribution of 
A av (BCS) = 74% of the average phenomenological gap , which is 1.42 MeV. 
The contribution of the terms containing the QPI terms to A a „ is equal to 
(-5)% of the average phenomenological gap . This result is similar to that 
for the ph channel in the sense that the contribution of terms corresponding 
to phonon exchange diagram and to the self- energy diagrams are opposite in 
sign, but in our case the contribution of the self-energy diagrams in the pp 
channel is small. 

Thus, we obtained that the total contribution of the QPI to the averaged 
phenomenological gap for 120 Sn is (31-5)=26% and the BCS-like part is 74%. 
One should emphasize , however, that these are just averaged figures and our 
method of determining the new pp interaction is the simplest one. But in any 
case we see that pairing in semi-magic nuclei is of a mixed nature - it is due 
to the BCS-like mechanism and the quasiparticle-phonon one. 

As it was mentioned above, for the calculations of nuclear physical charac- 
teristics with taking into account the QPI it is enough to know the new basis 
i\,A\. The question arises where one can see the effect of this new basis. 
The simplest reply is to compare the calculations with new and old basis. We 
have made this comparison (in the framework of our g 2 approximation) for the 
spectroscopipjiactors in 119 Sn and 121 Sn. The method of calculations has been 
described inull3. The results are given in Table 1. 

One can see easily that using our new (refined) basis improves mainly 
the description of the experiment. Of course, it would be desirable to find 
a better confirmation , i.e. a clearer manifestation of the effect . Probably 
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it can be seen in the calculations of other characteristics of low-lying levels 
of odd-mass nuclei. The experiments using the gamma-spectrometers of the 
EUROBALL-type can clarify this question. 



4 Strong coupling (g 2 > 1). Model calculations. 

Unfortunately, in this case it is necessary, generally speaking, to go beyond the 
Eliashberg-Migdal approximations, i.e. to add the vertex corrections to the 
rainbow ones. The reason is that, because our phonons are made up of the 
same qirasiparticles, there is no such a small parameter as in the Eliashberg 
theoryEl due to the existence of which these corrections are negligible. Actually, 
there can be some specific reasons for that in nuclei (see below) but , in any 
case, in finite nuclei we do not have such a small parameter g 2 as in the theory 
of ordinary superconductors. Thus, instead of the mass operators (^), we 
should take into account at least the following mass operators 

•••*•••••**•. 

M = n q + n— n o o (7) 

where the double line means the full Green function in its matrix form which 
includes G, G h , F^, F^l 

It is necessary now to solve Eqs. ^| with these mass operators. For the 
qualitative understanding and analysis of the case with strong coupling, here 
we consider a schematic 3-level model. We will take three levels with energies 
e=-l, 0, 1, each of them has the quantum number j — 11/2, the first one being 
occupied completely, parity of the low- and high-lying levels being positive and 
that of intermediate level negative. The occupation number of the intermediate 
level will be changed from to 2j + 1 = 12. In this model , to obtain A ^ 
0, the parameter of pp interaction should be G pp =0.035. The ph interaction 
was taken in our case to obtain the Bohr parameter a = g 2 /((2j + l)ui s ) (ui s - 
is the energy of phonon) equal to 2 for the most collective 2 + phonon. 

To estimate the contribution of the vertex correction we used the well- 
known results by Belyaev and Zelevinsky E3 according to which every vertex 

contains the factor w with a 6j-symbol u> = — (2j + 1) i ^. ^ I so that the 



J2j 

n-iteration contains a factor w n which decreases the contribution of the vertex 
correction as compared to the rainbow contribution. In our case w=0.76, so in 
what follows we restricted ouselves to w only. That means that at least the g 4 
terms are taken into account reasonably enough in our nonlinear model with 
the mass operators (^). For the beginning and simplification we have calculated 
the spectroscopic factors in our non-linear model with the mass operators (0) 
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Figure 1: Distribution of the single-particle strength calculated in the three-level model for 
the case with pairing, j=ll/2, a=2 ( see the text for details) The dispersion was taken as a 
smearing parameter for each of the calculations. 



, the above-mentioned parameter a being taken a=2, which corresponds to a 
realistic nuclear case. To obtain the distribution of the single-particle strength 
the iteration procedure was organized for the solution of the Dyson equation 
Eq. (|2|) with our mass operators. In our case this procedure must include also 
the refinement procedure for and Aj. 

The results are given in Fig.l. It turned out that it was necessary to do 
four iterations to obtain the convergency. This corresponds to the approxima- 
tion lquasiparticle® 4 phonons. We can see that the inclusion of the vertex 
correction is noticeable in the strength distribution. 

In particular, for the dominant level we obtained the decrease of the spec- 
troscopic factor from 0.49 to 0.40 due to the inclusion of the vertex corrections. 
The calculations have shown also that the role of the mass operators and 
M^ 2 ) for the levels far off the Fermi surface was small. 
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5 Conclusion 



We obtained that, if the simple procedures proposed are used to determine the 
new particle- particle interaction and to estimate the average effect according 
to Eq.(6), then the QPI contribution was 26% and the BCS-like mechanism 
gave 74% of the gap observed for 120 Sn. This means that at least for semi- 
magic nuclei pairing is of a mixed nature-it is due to the quasiparticlc-phonon 
and the BCS-like mechanisms, the first being mainly a surface one and the 
second mainly a volume mechanism. The effect of the QPI in the gap value 
can be observed probably in the experiments using the gamma-spectrometers 
of the EUROBALL-type to measure the characteristics of low-lying levels. 

For the case of nuclei with unfilled shells in both proton and neutron 
systems (strong coupling) our estimation has shown that the contribution of 
the vertex correction to spectroscopic factors is rather noticiable. One can 
think that this contribution is even more noticeable in transition probabilities 
for low- lying levels. 

S.K. thanks Prof. G.M. Eliashberg for useful discussions of the results. 
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